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Abstract 

By employing the perturbative QCD(pQCD) factorization approach, we calculated the partial 
next-to-leading order (NLO) contributions to Bg PP decays ( P = 7r,K,rj^'^ ), coming from 
the QCD vertex corrections, the quark-loops and the chromo-magnetic penguins, we found 
numerically that (a) for three measured decays Bg K~^7r~ , K^K~ and h~^'k~ , the consistency 
between the pQCD predictions and the measured values are improved effectively by the inclusion 
of the NLO contributions; (b) for Bg K^r]^'^ and K^tt^ decays, the NLO enhancements to the 
branching ratios can be significant, from ~ 50% to 170%, to be tested by the LHC experiments; 
(c) for the CP-violating asymmetries, the leading order pQCD predictions can also be changed 
significantly by the inclusion of the NLO contributions; (d) for Bg — > K^tt~ decay, the pQCD 
prediction for the direct CP asymmetry is Aqp{Bs — ^ K^tt~) = 0.26 it 0.06, which agrees very 
well with the only measured value available currently. 
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I. INTRODUCTION 



The two-body charmless hadronic decays of B or Bg meson are the good place to test 
the Standard Model (SM) and look for the signal of new physics beyond the SM. Since 
1999, more than 10^ events of BB pair production and decay have been collected and 
studied in the B factory experiments. In the Large Hadron Collider (LHC) experiments 
(ATLAS, CMS and LHC-b), besides those light Bu^ mesons, a huge number of heavier Bg 
meson production and decay events will be collected The study about the charmless 
decays of Bg meson is therefore becoming more interestirig then ever before. 

By employing the generalized factorization approachjl, 0] or the QCD factorization 
(QCDF) approach Q], about 40 Bg M2M3 (Mj stands for the light pseudo-scalar or 
vector mesons ) decay modes have been studied, for example, in the framework of SM 
[5I, 0| or in some new physics models beyond the SM |;7|. Many Bg meson decays have 
also been calculated, on the other hand, by employing the perturbative QCD (pQCD) 
factorization approach at leading order [sl, |9|, [lOl . 



Very recently, some next-to-leading order (NLO) contributions to some B — * M2M; 



decays have been calculated by employing the pQCD approach[ll|, Il2|, Il3|. One can see 
from those studies that the NLO contributions can change significantly the leading order 
(LO) pQCD predictions for some decay modes. It is therefore necessary to calculate the 
NLO contributions to those two-body charmless Bg meson decays, in order to improve 
the reliability of the theoretical predictions. 

we here focus on the calculations of NLO contributions to Bg PP decays ( P = 
71, K, T]^'^ ) in the pQCD approach. The NLO contributions considered here include: QCD 
vertex corrections, the quark-loops and the chromo-magnetic penguins. We expect that 
they are the major part of the full NLO contributions in pQCD approach [llj. The 
remaining NLO contributions in pQCD approach, such as those from the factorizable 
emission diagrams, hard-spectator and annihilation diagrams as illustrated in Figs. 5-7 in 



Ref. [13|, have not been calculated at present and should be studied as soon as possible. 

This paper is organized as follows. In Sec. [Tll we give a brief review about the pQCD 
factorization approach. In Sec. IIIIl we calculate analytically the relevant Feynman di- 
agrams and present the various decay amplitudes for the studied decay modes in the 
leading-order. In Sec. IIVI the NLO contributions from the vertex corrections, the quark 
loops and the chromo-magnetic penguin amplitudes are evaluated. We calculate and show 
the pQCD predictions for the branching ratios and CP violating asymmetries of Bg — > PP 
decays in Sec. V. The summary and some discussions are included in the final section. 



II. THEORETICAL FRAMEWORK 
A. Decay amplitude in pQCD 

In the pQCD approach, the decay amplitude is separated into soft ($AfJ, hard 
( H{ki,t) ), and harder( C{M]y) ) dynamics characterized by different energy scales 
i^QCD,t,mi,, Mw) The decay amplitude A{B — M2M3) can be written conceptu- 
ally as the convolution. 



2 



where kiS are momenta of light quarks included in each meson, and Tr denotes the trace 
over Dirac and color indices. C{t) is the Wilson coefficient evaluated at scale t. The 
hard function H{ki, k2, k2,,t) describes the four quark operator and the spectator quark 
connected by a hard gluon whose is in the order of AMb, and can be perturbatively 
calculated. The function $Af- is the wave function which describes hadronization of the 
quark and anti-quark in the meson Mj. While the hard kernel H depends on the processes 
considered, the wave function is independent of the specific processes. Using the wave 
functions determined from other well measured processes, one can make quantitative 
predictions here. 

Since the b quark inside the B meson is rather heavy, we consider the B meson at 
rest for simplicity. Using the light-cone coordinates, we define the emitted meson M2 
moving along the direction of n = (1, 0, Ot) and the recoiled meson M3 the direction of 
V = (0, 1, Ot). Here we also use Xi to denote the momentum fraction of anti-quark in each 
meson: 

Mb, , Mb , , Mb , 

Pb. = ^(1,1,0t), P2 = ^(1,0,Ot), P3 = ^(0,1,Ot), 

/Ci = (XiP+, 0, kiT), k2 = (X2P2+, 0, k2T), ks = (0, X3P3-, k3T). (2) 

Then, the integration over fcf , A;^, and k^ in eq.([T]) will lead to 



A ^ J dxidx2dx^hidhih2dh2h^dh^ 

•Tr [C{t)^B{xM^M,{x2M)^M-A^MH{x,M.t)St{x,)e-''^''^] , (3) 

where hi is the conjugate space coordinate of kiT- The large logarithms {\n.my/ /t) coming 
from QCD radiative corrections to four quark operators are included in the Wilson coeffi- 
cients C{t). The large double logarithms (In^ Xj) on the longitudinal direction are summed 
by the threshold resummation, and they lead to St{xi) which smears the end-point singu- 
larities on Xi. The last term, e~^^^\ is the Sudakov form factor which suppresses the soft 
dynamics effectively [13 ]. 



B. Effective Hamiltonian and Wilson coefficients 

For the studied Bg PP decays, the weak effective Hamiltonian He// for b 
transition can be written as lli 



^ q=u,c 1=3 J 

where Gp = 1.16639 x 10~^GeV~^ is the Fermi constant, and Vij is the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix element, Ci{n) are the Wilson coefficients evaluated 
at the renormalization scale fi and Oi{^) are the four-fermion operators. For the case of 
b —>■ d transition, simply makes a replacement of s by in Eq. (jlj) and in the expressions 
of Oi{ix) operators, which can be found easily for example in Ref. |15j. 

In PQCD approach, the energy scale "t" is chosen as the largest energy scale in the 
hard kernel H{xi, bi, t) of a given Feynman diagram, in order to suppress the higher order 
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corrections and improve the reliability of the perturbative calculation. Here, the scale "t" 
may be larger or smaller than the mh scale. In the range of t > rrib, the number of active 
quarks is Nf = 5, and the renormalization group (RG) running of the Wilson coefficients 
Ci{t) and LO and NLO level can be written as [l5|. 

= U{t,MwifC,{Mw)'''', 
Qitf'^'' = U{t,Mw,a),,C,{Mwf'^''. (5) 

The explicit expressions of C^^'^^°{Mw), the RG evolution matrix U{t,MwY°^ and 
U{t, M]Y, a) can be found easily, for example, in Refs. fl5| . 

In the range of fiQ < t < rrif,, the number of active quarks is Nf = 4, and we have 
similarly 

= ?7(t,m6,a),,C,(m,)^^o. (6) 

According to the analysis in Ref. we believe that it is reasonable to choose /io = 1-0 
GeV as the lower cut-off of the hard scale t, which is also close to the hard-collinear scale 
\/AmB ~ 1.3 GeV in SCET 16 1. In the numerical integrations we will fix the values 
Ci{t) at Cj(l.O) whenever the scale t runs below the scale /xq = 1-0 GeV. 



Wave functions 



As usual, we treat the B meson as a very good heavy-light system, and consider only 
the contribution of Lorentz structure 
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B, 



(7) 



with 



Ml x'^ 1 



where uji, is a free parameter and we take ujb = 0.5 ± 0.05 GeV for Bg meson. For a given 
Ub, the normalization factor Nb^ can be determined through the normalization condition 



fB 



(9) 



with Jb, = 230 MeV. 

For the r] — rj' system, we employ the quark-fiavor mixing scheme: the physical states 
7] and T]' are related to the flavor states rjg = {uu + dd) / \/2 and t]s = ss through a single 
mixing angle 0, 



cos (f) — sin I 
sin (h cos (h 



Vs 



(10) 
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The relation between the decay constants (/^, /^,, and {fq, fs, ), the chiral enhance- 
ment tUq and mg associated with the two-parton twist-3 rjg and 1]^ meson distribution 
amphtudes (DA's) have been defined in Ref.jllJ. The three relevant input parameters 
fq, fs and (p have been extracted from the data of the relevant exclusive processes |l7] : 



/, = (1.07 ± 0.02)/,, /, = (1.34 ± 0.06)A, 



39.3° ± 1.0°, 



with = 130 MeV. 

For the light pseudo-scalar mesons vr and^ as well as rjq and rjg, their wave functions 
are the same in form and can be defined as 18 



$(P,x,C) 
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=75 [f(p^{x) + mo0^(x) + Cmoif^f! - 1)0p(x)] 



(12) 



c 



where P and x are the momentum of the light meson and the momentum fraction of the 
quark (or anti-quark) inside the meson, respectively. When the momentum fraction of 
the quark (anti-quark) is set to be x, the parameter ( should be chosen as -f-l (—1). 
The expressions of the relevant DA's of the meson M = (vr, K, rjq, r]s) are the following 



^m{x) 
4'm{x) 



3/ 



M 



xil 



l + a'fcl'\t)+afCt'\t) 



2^6 . 
/m(1 



1 + ( 30r^3 - IpIi ] Cl'\t) 



2x) 



2^6 
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1 + 6 5r/3 - -ri^uj'i - -pL - tPm< (l - lOx + lOx^) 



2 ' 20' 5' 

{m^/m^,mK/m^,mqq/ml,mss/ml) for M = 
respectively [UJ, llS] . The Gegenbauer moments af^ have been chosen as 



(13) 
(14) 
,(15) 



with the mass ratio pm 



(tT, K, 7]q, r]s 



0.17, 



= 0, 



0.44, 



0.20. 



The values of other parameters are r/s = 0.015 and u 
polynomials C^it) in Eqs. (11311151) are defined as: 



(16) 

-3.0. At last the Gegenbauer 



Cf(t) = 3t, cr(t) = ^(3t^-l), C^,''{t) = ^{5e-1) 



3/2/ 



(17) 



1 /2 3 /2 

with t = 2x — 1. Like Ref. ^0|, we here also drop the terms proportional to C4 (t) in 
the LCDA's for (j)^, (j)^, and (j)^. 

There are many studies about the distribution amplitudes of light mesons and the 
relevant Gagenbauer moments (iSj. In recent years, the light-cone distribution amplitudes 
have been updated continually [l9|]. The inclusion of higher order terms and the variations 
of the Gegenbauer moments do affect the theoretical predictions for the branching ratios 
and CP-violating asymmetries of the considered decays, but the resultant changes of 



theoretical predictions are indeed not significant, according to the analysis in Ref. [10 
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(a) 



(b) 



(c) 



(d) 







(e) (f) (g) (h) 

FIG. 1: Feynman diagrams which may contribute at leading order to B'^ PP decays. 



III. DECAY AMPLITUDES AT LEADING ORDER 

The thirteen i?^ PP decays (P = vr, t], t]') have been studied previously in 
Ref. [lo[ by employing the pQCD factorization approach at leading order. The decay 
amplitudes as presented in Ref.[lo| are confirmed by our recalculation. In this paper, 
we focus on the calculations of some NLO contributions to these decays in the pQCD 
factorization approach. 

At the leading order, the relevant Feynman diagrams which may contribute to the 
— >■ PP decays are illustrated in Fig. [H For the sake of completeness, however, we 
firstly show the relevant LO decay amplitudes in this section based on our own analytical 
calculations. 

For 5° — > K^r]^'^ decays, the LO decay amplitudes are 



. K^r]) = ^(5° ^ K\) cos(0) - ^(5° 
K^) = ^(5° K\) sin(0) + ^(5° 



irV)cos(0), 



(18) 
(19) 



with 



-6 { fqFeK 2a3 + 04 - - -a-j + -ag - -a 



+MeK (C3 + 2C4 + 2Q 

+ fBsFaK 0,4 



1 



2 2 



10 



1 



c 
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+MaK (Cs - ^C,^ + M^k [c, - \C,^ } , 



(20) 
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+ (Me,. + M^k) (c, - iCg) + {M^l + Mj^k) [c, - \c,^ I , (21) 

where = 6 = Vt,V;^. 

For S° — > 7r°?7^') decays, the LO decay amphtudes are 

A{Bl ^ 7r%) = A{Bl ^ TT^g) cos(0) - A{B'l ^ 7r%,) sin(0), (22) 
A{Bl ^ 7r°V) = ^(4° ^ 7r°77g) sin(0) + ^(5° ^ 7r°r;,) cos(0), (23) 

with 

A{Bl ^ 7r°r/,) = {fsF^r,^ a, + M„,„ C^) 

-\i't [fB,Fa,„ (ay + ag) + M„,„ Cio + Cs] , (24) 
^ 7r\) = (/.i^e.. a2 + M,,, C2) 

-^e; [/.i^e,. (ag - ar) + M„,. (Q + Cio)] , (25) 

where ^ = K^KI,. CI = Vt,V*. 

For S° ^ r]rj, rjr]', rjrj decays, the LO decay amphtudes are 



^PiA{Bl - 




= cos^ (0)^(5° ^ 


77,77,) - sin(20)^(5° ^ 


VqVs) 








+ sm'{cl>)AiBl 


?7s?7s), 




(26) 


A{Bl - 


rjr]') 


= [^(4^^77,77,) 


- ^ 77,77,)] C0S((/)) 


sin(0) 








+ cos(2(/))^(S° 


^ ?7g?7s), 




(27) 


V2A{Bl - 




= sin^(0)^(5° ^ 


77,77,) +sin(20)^(S°^ 


VqVs) 








+ cos' (0)^(5° 


VsVs)- 




(28) 
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with 



(29) 



1 



fnFeris I — 2a5 — -ay + -ag 



f 111 

fsFerj, ( as + 04 - as + -ay - -ag - -aio 



(30) 



+ (Me,, + M„,J (^Ca + C4 + Ce - ^Cg - ^Cg - 

For 5° — i> i^T+TT" , i^r°7r° K^K~ and K^K^ decays, the LO decay amphtudes are 
A{K+t:-) = e« (/.ri^ei^ oi + Mei^ Ci) - {/.FeK (04 + Oio) + UF^^ (ae + ag) 

+MeX (C3 + C9) + fBsFaK {^A " ^^lo^ + - ^Cg 



(31) 



(32) 



V2^(X°7r°) = e« {UFeKa2 + M^i^Cs) - 6 I^UFex (^-ai - ^ay + ^ag + ^aio 

- (A^e^^ + fsF^^) (a, - ^ag) + Me, (^-^3 + ^Cg + ^Cg + ^Cio) 
(^a, - ^aio) - M„fe - ^Cg^ - M^k ([^5 - ^^^7) } , (33) 



A{K+K-) = C{fkFekai + M,kCi + MakC2) 

-C't {fkFek (04 + oio) + fkF^ (oe + as) 

+Mefe (Cs + Cg) + M,^^ {C, + Cy) + /^^Fj^^ (^ag - ^ 



ag 



+M„fe ( C3 + C4 - -Cg 



;Cin 1 + M^^ ( C5 - ^Cy 



1 

-( 
2 



<^ (^6 - ICs^ + [Mak (Q + C,o) + m5 (Co + Cg)] | , (34) 
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+Mak C3 + C4 - -Ca - -a 
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10 



+ 



m2 ( c, - -Cs 



(35) 



The —>■ TT+TT and —>■ tt^tt^ decays are pure annihilation processes, and the LO 
decay amphtudes can be written as: 



■ (36) 



The individual decay amplitudes, such as F^k, Feri^,etc., appeared in Eqs. fll8ti3ip . are 
obtained by evaluating the corresponding Feynman diagrams analytically. For the Bg —>■ 
M3 transitions, where the meson M3 absorbed the spectator s quark, the corresponding 
decay amplitudes can be written as 



SnCpM^^^ / dxidxs / 61^6163^63 61) 
° Jo Jo 

■ { [(1 + x,)<Piixs) + r3(l - 20:3) (0f (0:3) + (Plixs))] 

■as(tl)he{xi,X3, 61, 63) exp[-Sab(tl)] 
+2r30f(x3) ■ asitl)K{xs,Xi,b3,bi) exp[-Sabitl)]} , 



(37) 



^eMs = lQTrCFM^j2 / dxidxs / bidbib3db3 4>Bs{xi,bi) 

Jo Jo 

■ { [0^(X3) + r3(2 + X3)0f (X3) - r3X30f (X3)] 

■asitl)heixi,X3, 61, 63) exp[-Sabitl)] 
+2r30f(x3)as(tg)/ie(x3,xi, 63,61) exp[-S'a6(te)]} , (38) 



M, 



32 



dxidx2dxs / 61(^6162^62 0B^(xi, 61)0^(^2) 



■ { [(1 - X2)(l>t{Xs) - r3X3 (0f (X3) - 01 (X3))] 
■as(^e)^n(2^1> 1 - ^2, X3, 61, 62) exp[-Scditl)] 

- [{x2 + X3)03 (X3) - r3a;3 (0f (^3) + ^^(a^s))] 
■as{4)hn{xi, X2, xs, 61, 62) exp[-Scd{tt)]} , 



(39) 
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bidbib2db2 (f)B,{xi,bi) 



•{[(1-X2)0^(X3) (0r(a;2)+0^(a;2)) 

+r3X3 {4)2 ix2) - (t)lix2)) (0f (2:3) + 03 (2:3)) 
+r3(l - X2) (0f (X2) + 0f(.T2)) (0f (X3) - 01(^3))] 

■as{tl)hn{xi, 1 - a;2, X3, 61, 62) exp[-S'cd(te)] 

+ [X203 (2:3) (0f (a;2) - 0n^2)) + ^32^2 (^^SJa) - 02 (^2)) (^f (^^s) " ^H^s)) 

+r3a;3 (0^(0:2) + 0^(2:2)) (0^(0:3) + 0r(x3))] 
■as{tl)hn{xi,X2,X3,bi,b2)eiq)[-Scd{tt)]} ■ (40) 

32 

^eMs = —;='^CfM^^ / dxidx2dxz / &irf&i&2'^&2 0B,(a^i, &i)0^(a^2) 



aMg = SttC^M^^ / dx2dx3 / b2db2b:idb:i • { [(2:3 - l)0^(a;2)03 (xs) 




• { [(a;2 - X3 - 1)03 (2^3) + r3a;3 (0f (^3) + 01(0:3))] 



■as(te)/?.n(.«i, 1 - X2, .T3, 61, ^2) exp [-S'crf(ie)] 
+ [a;20^(a;3) + r^x^ (0f (1X3) - 01(2:3))] 



•Qis(^e)^n(^l' ^2, 2^3, ^2) exp[-S'cd(te)] } , 



(41) 




-4r2r30^(a;2)0^(a:3) + 2r2r3a:30^(a;2) (0f (^3) - 03 (2^3))] 

■as{tl)ha{x2, 1 - 2:3, 62, ^3) exp[-S'e/(te)] 
+ [a;20^(a:2)0^(x3) + 2r2r3 (0f (^2) - 0^(x2)) 0f (xg) 
+2r2r3a;2 (0^(x2) + 0^(0:2)) 0f (0:3)] 



(42) 




• { [2r20f (0:2)0^(0:3) + (1 - X3)r30^(x2) (03^(0:3) + 0^X3))] 

■as{tl)ha{x2, 1 - 0:3, 62, ^s) exp[-5'e/(te)] 
+ [2r30^(x2)0f (xs) +r2a:2 (0f (2:2) - 0^(2:2)) 03 (2^3)] 

■as{tl)ha{l - X3,X2MM) Q'^V[-Sef{tl)\] , 



(43) 
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32 f°° 
MaMg = —^ttCfM^^ / dxidx2dx3 / bidbib2db2(j)Bsixi,h) 
v6 Jo Jo 

■ { [-a;20^(a;2)03 (^3) - 4r2r30f (xa)^^ (xg) 
+r2r3(l - X2) (0f (X2) + 0na;2)) (0^(2:3) - (Plix-i)) 
+r2r3X3 (0^X2) - (f)2{x2)) (0f (3:3) + ^n^^s))] 
■as(4)^na(a;i, X2, X3, bi, 62) exp[-5'g/,(tg)] 

+ [(1 - X3)0^(X2)0^(X3) + (1 - X3)r2r3 (0^(0:2) + 4>^{X2)) (^f (X3) - (Xs)) 

+X2r2r3 (0^X2) - 0^(x2)) (0f (0:3) + 0^(3^3))] 
■"^(^e)^na(a;ba;2,a;3,&i,&2)exp[-5'<,/,(tg)]} , (44) 



/■I /-oo 



= —nCpM^^ / dxidx2dx3 / 6io?6i62C?&2 &i) 
v6 Jo ^0 

■ { [r2(2 - X2) (0r(^2) + 0^(X2)) 0^(X3) + r3(l + X3)0^(X2) (0f (X3) - <PI{X3))] 
■asitl) exp[-Sgh{tl)]hna{Xi, X2, X3, bi, 62) 

[r2X2 (0f (X2) + 4>^{X2)) (f)i{x3) - (1 - a;3)r30^(x2) (0^(0:3) - ^3 (2^3))] 
s(^e)^na(a;i5a;2,a;3,&i,&2)exp[-5g/,(t^)]} , (45) 



+ 

•a. 



= —ttCfM^ / dxidx2dx3 / bidbib2db2(t)Bs{xi,bi) 
v6 "Jo io 

• { [(X3 - 1)0^(X2)0^(X3) - 4r2r30f (X2)03''(a^3) 

+r2r3X3 (0f (X2) + 02 (a;2)) (0^(3:3) - 03 (X3)) 
+r2r3(l - X2) (0f (X2) - 0^x2)) (0f (0:3) + 0r(x3))] 

■as(tl)hna{xi, X2, X3, 61, 62) exp[-S'g/i(tg)] 
+ [x20^(x2)03 (X3) + X2r2r3 (0^ (^2) + 02 (2:2)) (0r(a:;3) - (I>l{x3)) 

+r2r3(l - X3) (0f (X2) - 0^(X2)) (0f (X3) + (P^{X3))] 

■(^s{tl)Ka{xi,X2,X3,bi,b2)exp[~Sgh{tl)]} , (46) 
where Cp = 4/3 is the color-factor, 

r2 = K^VMb^ and r3 = VTi^'^ 1Mb, with the chiral 
enhancement factor mo for meson M2 and M3. Here (FeMa, -f^ilfg) and (MeMs, -f^eMa ) come 
from the factorizable emission diagrams ( Fig. la and lb) and the non-factorizable hard 
spectator diagrams (Fig.lc and Id), respectively; and {FaMz-,PaM-^ ^^^1 [MaM^iPalh'^^ 
are obtained by evaluating the factorizable annihilation diagrams ( Fig.le and If) and 
the non-factorizable annihilation diagrams (Fig.lg and Ih), respectively. The explicit 
expressions of the hard energy scale (tl,tl, ■ ■ ■ , tg), the hard functions {he, hn, ha, h^a, h'^^), 
the Sudakov factors {Sah{te), Scd{te), Sef{te), Sghite)) can be found in Appendix A. 

IV. NEXT-TO-LEADING ORDER CONTRIBUTIONS 

At next-to-leading order, firstly, the NLO Wilson coefficients Ci{Mw), the NLO RG 
evolution matrix U{t,m,a) |15|], and the as{t) at two- loop level will be employed. 
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ivL X : X M 




(b) 



(f) 



(c) 




(g) 



(d) 



FIG. 2: Feynman diagrams for NLO contributions: the vertex corrections (a-d); the quark- loop 
(e-f) and the chromo-magnetic penguin contributions (g-h). 

Secondly, the NLO hard kernel H^^\a'^) should be included. All the Feynman dia- 
grams, in other words, which lead to the decay amplitudes proportional to a'^it), should 
be considered. Such Feynman diagrams can be grouped into following classes: 

I: The vertex corrections, as illustrated in Figs. [2^-l2li, the same set as in the QCDF 
approach. 

II: The NLO contributions from quark-loops, as illustrated in Figs. [2^-[2f. 

Ill: The NLO contributions from chromo-magnetic penguins, i.e. the operator Ogg, 
as illustrated in Figs. [2^-l2Ji- There are totally nine relevant Feynman diagrams as 
given m Ref. 0, 

if the Feynman diagrams involving three-gluon vertex are also 
included. We here show the first two only, and they provide the dominant NLO 
contributions, according to Ref. [20]. 

IV: The NLO corrections to the LO emission diagrams(la,lb), the LO hard- 
spectator(lc,ld) and the LO annihilation diagrams (le-lh), as illustrated in Fig. [31 
One can draw more than one hundred such Feynman diagrams in total, but we here 
show representative ones only. 

At present, the calculations for the vertex corrections, the quark- loops and chromo- 
magnetic penguins have been available and will be considered here. We expect that they 
are the major part of the full NLO contributions in pQCD approach For the Feynman 
diagrams as shown in Fig. 3, however, the analytical calculations have not been completed 
yet. Of course, these Feynman diagrams should be calculated as soon as possible, in order 
to improve the reliability of the pQCD predictions. 



A. Vertex corrections, quark-loops and chromo-magnetic penguins 

The vertex corrections to the factorizable emission diagrams, as illustrated by Figs. 2a- 
2d, have been calculated years ago in the QCD factorization approach^, 0]. For Bs — > 
PP decays, the vertex corrections can be calculated without considering the transverse 
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FIG. 3: The typical Feynman diagrams, which provide NLO contributions to the factorizable 
emission amphtudes (al-a4), the hard-spectator amphtude (bl-b4), and the annihilation ampli- 
tudes (cl-c4) 



momentum effects of the quark at the end-point [11], one can use the vertex corrections 
as given in Ref. 0] directly. The vertex corrections can then be absorbed into the re- 
definition of the Wilson coefficients ai{fi) by adding a vertex-function Vi{M) to them 



ai(/i) 
aj(/i) 



+ ^^Cp^Vm. for z = 1, 2; 



aj(/i) + 



An 



-C 



An 



-C. 



VjiM), for J = 3, 5, 7, 9, 
V,{M), for J =4, 6, 8, 10, 



(47) 



where M is the meson emitted from the weak vertex. When M is a pseudo-scalar meson, 
the vertex functions Vi{M) are given ( in the NDR scheme) in Refs. [g], 11|: 



Vi{M) 



121n^-18 + ^/;rfx0^(x)(7(x), for « = 1-4,9,10, 
421n^ + 6-^/,'dx(/.^(a;)(7(l-x), for ^ = 5, 7, (48) 
-6 + ^ /o dx(l)^j{x)h{x), for i = 6, 8, 



where /m is the decay constant of the meson M; (j)M{x) and (pM^^) the twist-2 and 
twist-3 distribution amplitude of the meson M, respectively. The hard-scattering functions 
g{x) and h{x) in Eq. fH5]) are: 



h{x) 



3 mx — tn 



1 



X 



+ 



2 Inx 



2Li2[x) — In X + 
2Lz2(x) - In^ X - (1 + 2in) Inx 



(3 + 2in) Inx — (x ^ 1 — x) 
X ^ 1 — x), 



(49) 
(50) 



13 



where Li2{x) is the dilogarithm function. As shown in Ref. [llj, the /x-dependence of 
the Wilson coefficients ai{fi) will be improved generally by the inclusion of the vertex 
corrections. 

The contribution from the so-called "quark-loops" is a kind of penguin correction with 
the four quark operators insertion, as illustrated by Fig. [2^ and[2]f. We here include quark- 
loop amplitude from the operators Oi,2 and O^-q only. The quark loops from Oy-io will 
be neglected due to their smallness. 

For the b ^ s transition, the effective Hamiltonian Hfjj- which describes the contribu- 



tions from the quark loops can be written as |11 



rriql) 
^eff 



-EE 

q=u,c,t ql 



Gp 

7! 



Vq,v;, 



(.-7p (1 - 75) T%) (qYT-q') , (51) 



where P is the invariant mass of the gluon, as illustrated by Figl2fe. The functions 
C^''\lJ,,P) are given by 



for q = u,c and 



GW(/x,/2)- 



C2(/i) 



(52) 



CM+ G'-'''\fi,l')[C,{fi) + CM] (53) 



q' =u,d,s,c 



for q = t. The function G^''\jj,,l'^) for the loop of the q{q = u,d,s,c) quark is given by 

Q 



G(«)(/i,/^ 



-4 / dxx{l — x) In 
Jo 



— x(l — x)f 



(54) 



m. 



is the possible quark mass. The explicit expressions of the function G^i\fx,f) after 



the integration can be found, for example, in Ref. [U 



The magnetic penguin is another kind penguin correction induced by the insertion 
of the operator Osg, as illustrated by Fig.2g and 2h. The corresponding weak effective 
Hamiltonian contains the b ^ sg transition can be written as 



Tj-mp 



with the magnetic penguin operator 

9 



G 



V2 



m, d, a^^ (1 + 75) ^ Gl, b, 



(55) 



(56) 



where i, j being the color indices of quarks. The corresponding effective Wilson coefficient 
Gll^ = Gsg + C5 [11 • 



It is worth noting that there are totally nine Feynman diagrams involving the Ogg 
operator as shown in Ref. 2^, if the Feynman diagrams involving three-gluon vertex are 
also included. We here show the first two only, say Fig.2g and 2h, and they provide the 
dominant NLO contributions from the O^g operator, according to Ref. [20| . 
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B. NLO decay amplitudes 

Now we are ready to calculate the decay amplitude corresponding to the quark-loops 
and magnetic penguins. By analytical evaluations, we find two kinds of topological decay 
amphtudes A^M2M3 •^M2M3' respectively. For B — > rjnK'^ decay, we find 

^vIko = --^m%Cl j dxidx2dx3 j bidbib3db3(f)BSxi) {[{1 + X3)(p^Jx2)(t>iix3) 

+2r,,„0jl(a;2)0i(a;3) + rxil - 2x3)0,,Ja;2) [(pxixa) + 4>k{x3)] 

+'^VK(P^^{X2) [{2 + X3)(j)^{x3) -X30]^(X3)]] 

■al{tl)K{xi,xs, 6i, 63) exp[-Sai,{tl)]C^''\tl f) 

l{tl')Kix3, xi, 63, 61) exp[-5„,(4°)]C('') (^^^ I")} , (57) 



M^f^o = ^ ^ dx,dX2dX3 b^db^b2db2b3db3(t>BX^,){[-{l-X3)[(t>i{x3) 

+rK{^(p^{x3) + 0]^(x3)) + ri^a;3 (0^(0:3) - 0^(a:3))]0rt (2^2) ~ ^'(n 2^2(1 + ^3) 
x(30j'„(a;2) -0j:(a;2))0^(x3) -r,„r^(l -a;3)(30jl(x2) + 0j:(x2))(0^(a;3) 
-0x(a:3)) - VKX2{1 - 2a;3)(30jl(x2) - <„(:r2)) (0^(0:3) + ^^(xa))] 
■a,2(t^)/i^(xi, ,X2, ,X3, 61, 62, &3) exp[-5„KtD]C|f (^e) 
- [ATKct>i{x2)ct>'K{^3) + 2r,„r^a;2(3<(x2) - 0^„(x2)) 0^(0:3)] 

•a^(i^°)/i;(xi, X2, X3, 61, 62, 63) exp[-5e<i(i^°)]<^(if ) } ■ (58) 



A^S^^o,, = ^|-^ y dXidX2dXs J b,db,bsdbscj>BXx,){[{l + Xs)cPi{x2)cj>^Xxs) 

+2rK(l>K{x2)(l>i{x3) + r^,(l - 2x3)0i^(x2) K(^3) + 0j.(a;3)] 
+2r,^r^, 0^(2:2) [(2 + 0:3)0^ (0:3) - X30^^(x3)]] 
■alitl)Kix^, X3, 61, 63) exp[-Sab{tl)]C^^\tl Z^) 
+ [%>^(^2)0j(a;3) +4r;^r^>^(x2)0j(x3)] 

•a,^(t^°)/ie(x3, x„ 63, &i) exp[-5„,(tf )]C('^)(t^°, r^)} , (59) 
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16m 



f 




■ { [-{1 - xs) [2<(X3) + r,, [3<(a;3) + ^(xs)] 
+V^{(l>i{xs) - (j^lixs))] 0^(0:2) 
-rKX2{l + X3) [3(j)^{x2) - 0^^(2:2)] (p^S^s) 
-rKr^,(l - X3){3(j)^{x2) + 4>k{x2)){(P^S^3) 

-^(•^3)) ~ ^KV^C^ - 2X3){3(P^{X2) - 4'K{x2)){(pi{x3) + ^Jll^^s))] 

■a''M)Cli^ {tl)hg{x,, X2, X3, 61, 62, bs) eM-Scditl)] 

■al{tf%{x,, X2, xs, 61, 62, h) exp[-5e,(4°)]C|f } , (60) 



= -16m%^-^ / dXidX2dX3 / 6id&i63^i&3 0B.(xi){[(l+X3X(^2X(x3) 



+2r^Jl{x2)<j)l{xs) + r,,(l - 2x3)0,. (0:2) (0^(0:3) + (J^lixs)) 
+2rr,JnA^^{x2){{2 + ,X3)0j(a;3) - x^c^J^^ix^))] 
■al{te'^)K{x,, Xs, h, 63) exp[-Sab{tl)]C^'\tl, f) 
+ [2r,>J(a^2)0j(a:3) +4r,^r^>j;(a:2)0j(a:3)] 



M^]^^ = -32m|^-^ / dx,dx2dx3 / 6id6i62(^&2M&3 0B.(a:i) { [(1 - 2:3) [20;^^(a;3) 



(30^(0:3) + 0jl(a;3)) + V3i(l)iix3) - 0jl(2;3))] </';t(^2) 

-r^.X2(l + a;3)(30,^^(a;2) - (p'^S^2))(j)^S^3) 

- XS){3<PI{X2) + 0J,(X2))(0J(X3) - ^liXs)) 
-rr^sV2{l - 2X3)(30J^(X2) - 0J^^(X2))(0J(S3) + 0jl(2;3))] 

•a^(i^)/i,(xi, X2, xs, h, b2, bs) eM-Scd{tl)]Cli^itl) 

+ i^^VsKi^^^Ki^s) + 2rr,Jr,sX2mi{x2) - 0j^(x2))0j'^ (Xg)] 

■al(tl')h'^(x,,X2, xs, 61, 62, bs) eM-Scd{tl')]Clf/{tl') } . (62) 



It is easy to see that the decay modes -B^ T^^Tjn, 'n'^Tjsjrjn'rin and ?7„??s, do not receive 
the NLO contributions from the quark-loop and the magnetic-penguin diagrams. 

For other Bg PP decays involving no r]^'^ mesons, the corresponding quark-loop and 




■al{t'^)K{xs, x„ bs, 61) exp[-Sa,{tl')]C(''\tl\ l")} , 



(61) 
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magnetic penguin amplitudes are of the form 

•^!r-i^+ = ~^^B,^J= J dXidX2dX3 J bidbib3dh(l)B,{xi){[{l+X3)(j)^{x2)(l)K{x3) 
+2r^0^(x2)0^(x3) + rx(l - 2x3)0^(x2)(0^(x3) + (plixs)) 

+2vk(I)^{x2){{2 + X3)(f)^ix3) - X3(f)l{x3))] 

■al{tl)K{x,, xs, 61, 63) eM-Sab{tl)]C^''\tl, f) 
+ [2rK(l)^{x2)(l)K{x3) + 4r^rif0^(x2)0^(x3)] 

■a'M')K(x3, x^, 63, 61) eM-Sa,(tl')]C^'^\tl', l")} , (63) 



= -16m^^^^ / dxidx2dx3 / bidbib2db2b3db3(j)B,{xi) 
V ^-'Vc ^0 ^0 

• { [(1 - X3) [2<f>i{x3) + rx(30^(x3) + 0^(X3)) 

+rKX3{(l)K{x3) - (f)Kix3))] (f^ni^^) 

-r^X2{l + X3){3(j)^{x2) - (j)l{x2))<j)K{x3) 

-Vk(1 - X3)(30^(X2) + 0^(X2))(0:^(X3) - 0^(2^3)) 

-VkX2{1 - 2X3)(30^(X2) - 0^(X2))(0^(X3) + 0^(X3))] 

•Q;^(i^)/lp(xi, X2, X3, bi, 62, &3) exp[-5e<i(^e)]C8fl'^(^e'') 

+ [4rK(l>^{x2)(l>^ix3) + 2r^rKX2iS(/>^ix2) - 0^(X2))0^(X3)] 

a'M>g{^u X2, X3, 61, 62, 63) exp[-5,,(t^°)]C|f (4°) } , (64) 



= (66) 



C ^ /"I /"OO 

= -8m^^— ^= / dxidx2dx3 / 6i(i6i63(i63 ^^^(xi) 

v^J-^cJo Jo 

• { [(1 + a;3)0^(a;2)0^(a;3) + 2rfc(/.f (a;2)0i(a;3) 

+rK(l - 2a;3)0fc(a;2) (0^(0:3) + 0^(2:3)) 

+2rfcrK0r(a;2)((2 + a;3)0^(x3) - X34>Kix3))] 

■al{tl)K{xu X3, 61, 63) exp[-Sab{tl)P'\tl, 

+ [2r;^0^(x2)0^(x3) + ArkrK(t>l{x2)^K{x3)\ 

■a'M>e{x3, x„ 63, 61) exp[-5„,(t^°)]C('^)(t^°, T^) } , (67) 
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= -lQm%^^^ / dxidx2dx3 / hidhih2dh2hdh^(t)BX^i) 
v^J-^cJo Jo 

■ { [(1 - 0:3) [20^(0:3) + ri,(30^(x3) + <Pl{x3)) 

-rKX2{l + a;3)(30^(x2) - 0]^(a;2))0^(a;3) 

-4(1 - X3)(30^(X2) + 0^(a;2))(0^(a;3) - 0^(a;3)) 

-4x2(1 - 2x3)(30^(x2) - 0^(a;2))(0^(x3) + fA^z))] 

■a'M)hg{xi, X2, X3, 61, 62, 63) exp[-5,rf(t:!)]C|f (t^) 

+ [4rii-0^(x2)0^(x3) + 2r/^r/^X2(30^(x2) - 0K(a;2))</>]^(a;3)] 

■a'M'Hi^u X2, X3, &i, &2, h) exp[-S,,{tl')]Cf/itl') } , (68) 



Mfo,o = Mit., (69) 
•^1^.0 = M[%^. (70) 

The functions he,hg,h'g, the hard scales tg and tg°, the Sudakov factors Sab{t) and Scd(t), 
appeared in Eqs. (|Fr|) - (|^ . will be given in Appendix A. 

For the pure annihilation-type decays — > tt^tt" and vr^vr^, they do not receive 

the NLO contributions from the vertex corrections, the quark loops and the magnetic- 
penguins. For Bg —>■ 71^1] and B^ tt^t]' decays, the quark-loops and magnetic penguins 
also do not contribute. For other decays, the NLO contributions will be included in a 
simple way: 



q=u,c,t 
q=u,c,t 
q=u,c,t 

A^-K+ ^ Yl ^M'~K++^tMl'l^^, 

q=u,c,t 
q=u,c,t 

Ak-k+ Ak-k+ + Y, + 

q=u,c,t 

AkOko AkOko + ^ ^qM^'llj^o + ^'tM^llj^o, (71) 



q=u,c,t 

where = VgbV*^, and = VgbV*^ with q = u,c, t. 
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V. NUMERICAL RESULTS 



In the numerical calculations the following input parameters will be used. 

A^^ = 0.225GeV, /b, = (0.23 ± 0.02)GeV, fx = 0.16GeV, 

= O.lSGeV, Mb, = 5.37GeV, ttik = 0.494GeV, = 547.9MeV, 
= 0.958GeV, = 1.4GeV, = 1.9GeV, tbo = 1.470 ps, 

TTib = 4.8GeV, Mw = 80.42GeV. (72) 

For the CKM matrix elements, we also take the same values as being used in Ref . [13] , 
and neglect the small errors on Vud, Vus, Vts and Vtb 

iKdl = 0.974, iKsl = 0.226, = (3.68;°;^^) x 10-^ 

\Vtd\ = (8.20;°:^?) X 10-3, \Vts\ = 40.96 x 10~^ \Vtb\ = 1.0, 

a={99tLr, 7 = (59.0l|^r, arg [-1^,,^;,] = 1° . (73) 



A. Branching Ratios 



For the considered PP decays, the decay amplitude for a given decay mode with 

h d transition can be generally written as 



^(s° /) = v^bv:,T - Vtbv:,p = Vutv:,T [i + 



ze 



i{-a+5) 



(74) 



where a is the weak phase (one of the three CKM angles), 6 = arg [P/T] is the relative 
strong phase between the "T" and "P" part, and the parameter "z" is defined as 



(75) 



The ratio z and the strong phase 6 can be calculated in the pQCD approach. The CP- 
averaged branching ratio, consequently, can be defined as 



Bt(B 







/) 



32'KmB 



(76) 



where r^o is the lifetime of the Bg meson. 

For the case of 6 — > s transition, we have similarly 

^(4° ^ /) = V^bV:j' [1 + z'e'^^+'^] 



(77) 



where 7 is also one of the three CKM angles, 6 = arg[P'/T'] is the relative strong phase, 
while the parameter "2"' is defined as with 



P' 



T 



(78) 



In Table [T] we show the pQCD predictions for the CP-averaged branching ratios of the 
thirteen Bs —>■ PP decays. The label LO means the leading order pQCD predictions. The 
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TABLE I: Branching ratios (xlO~^) of PP decays in the pQCD approach. The label LO 

means the leading order pQCD predictions, while +VC, +QL, +MP, as well as NLO means 
that the vertex corrections, the quark loops, the magnetic penguins, and all the above NLO 
corrections are added to the LO results, respectively. The errors in the table are defined in the 
context. For comparison, we also cite the leading-order pQCD predictions as given in Ref. 
the QCDF results in Ref. jsl, and currently available data 



Mode 


Class 


LO + VC + QL + MP 


NLO 


pQCD[10] QCDF[6] 


Data 






C 


0.10 


0.08 


0.12 


0.10 


n 1Q+0.03+0.02+0.02 
"•^^^-0.06-0.00-0.04 


n 1 1+0.08 

"•-LJ--0.il 


'-'•"J^-0.33 




B^ ^ 




C 


0.70 


1.09 


1.13 


1.29 


J^-O'-0.51-0.11-0.21 


rv 79+0.36 
"• ' ^-0.24 










Pew 


0.04 


0.03 






r, r,Q+0.01+0.00+0.00 

"•'^■-•-0.00-0.00-0.01 


n5+°-°2 


0.081°:°^ 








Pew 


0.09 


0.08 






n ns+0 03+0 00+0 01 

"•"°-0. 02-0. 00-0. 01 










rjrj 


P 


7.4 


7.8 


9.1 


12.0 


in n+3-3+0.0+0.6 

^"■"-2.5-0.03-0.8 


8.0l|t 


15.6^ 






r]r]' 


P 


21.6 


31.7 


28.6 


37.4 


OA Q+11. 3+0.03+2.8 
"^^•^-8.6-0. 1-4.1 


21.0+1^ 


54 n+^2.8 

'^^•"-29.1 




B^ ^ 


rj'r]' 


P 


14.9 


30.3 


18.9 


27.1 


oc: 9+8.1+0.1+1.9 
^^•^-6.1-0.0-2.3 


14.0^1:? 


41 7+47.5 
^-■-^ ' -24.9 




B^ 


K+TT- 


T 


7.0 


6.3 


6.1 


6.0 


n O+2.5+0. 5+0.4 
"•"J-1. 8-0.5-0.3 


7.6+^1 


10.2160 


5.0 ± 1.3 






C 


0.16 


0.30 


0.18 


0.18 


n 9c;+0.09+0.03+0.02 
"•^'-'-0.06-0.01-0.03 


16+°- ^2 

"•-'-'-'-0.07 


n 49+0.62 

"•^^-0.35 








ann 


0.70 








A C7+0.14+0. 01+0.20 
"•'-"-0.11-0.00-0.19 


"•'^'-0.16 


"■"^-0.02 


0.53 ±0.5 


B^ ^ 




ann 


0.35 








n 90+0.07+0.01+0.10 
"•^^-0.06-0.00-0.10 


n 98+0.09 


"•"-■--0.01 




B^ 


K+K- 


P 


11.8 


15.5 


16.0 


15.4 


1 ^ R+S. 0+0.7+0. 8 
-^'J-"-3. 8-0.3-0. 7 


13.6l|6 


99 7+27.8 
' -13.0 


24.4 ±4.8 


B^ 




P 


14.3 


17.2 


18.0 


17.5 


IS n+4-6+0.0+0.7 
-^"•"-S. 9-0.0-0. 6 


-ic: f;+9.7 
-'-'^•"-6.0 


24 7+2^-^ 

' -14.0 





label ±VC, ±QL, ±MP, and NLO means that the vertex corrections, the quark loops, 
the magnetic penguins, and all the considered NLO corrections are included, respectively. 
The errors as shown for the NLO pQCD predictions correspond to the uncertainties of the 
various input parameters. The first error comes from ujb = 0.50±0.05 and fs^ = 0.23±0.02 
GeV. The second error arises from the uncertainties of the CKM matrix elements \Vub\ 
and \ Vcb\, as well as the CKM angles a and 7 as given in Eq. (175]) . The first two errors are 



defined in a similar way as that in Ref. |10 



The third error comes from the uncertainties of 

+0.10 



relevant Gegenbauer moments: af = 0.17±0.05, af = 0.20±0.06 and = 0.44+^;^^. We 
here assign roughly a 30% uncertainty for Gegenbauer moments to estimate the resultant 
effects on the theoretical predictions of the branching ratios. 

For the sake of comparison, we also list the leading order pQCD predictions as given 
in Ref. [lo| and the theoretical predictions based on the QCD factorization approach 
in Table [11 The corresponding errors of the previous LO pQCD and QCD F predictions 
denote the combined error: the individual errors as given in Refs. Ja llOf are added in 
quadrature. The currently available experimental measurements |2ll . |22| . |23| are also 
shown in the last column of Table [B 

From the numerical results about the branching ratios, one can see that 

• The LO pQCD predictions for branching ratios of Bg —>■ PP decays as given in 
Ref. M are confirmed by our independent calculation. The very small differences 
are induced by the different choices of the scales Aq^^ and Ag^^: we take A^^^ = 
0.225 GeV and the corresponding Ag^^, = 0.287 GeV, instead of the values of 
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AgJ.^ = 0.193 GeV and aJ^J.^ = 0.25 GeV as being used in Ref. [10 



In this paper, the NLO contributions are taken into account partially. The consid- 
ered NLO contributions can interfere with the LO part constructively or destruc- 
tively for different decay modes. For most decays the changes of the LO results 
are moderate and reasonable. The theoretical uncertainty from uj^ = 0.50 ± 0.05 is 
dominant, while the error from the uncertainty of CKM elements is small. And the 
total theoretical error is in general around 30% to 50% in size. 

For the "tree" dominated decay Bs K^tt^, the NLO pQCD prediction agrees 
with the data within one standard deviation. The agreement between the pQCD 
prediction and the measured value is improved due to the inclusion of the considered 
NLO contribution. 

For the three "Color-suppressed" decays, the NLO enhancement can be significant, 
from ~ 50% for the K^r] and K^t^^ decays to ~ 170% for the Bs K^r]' 

decay. The differences between the LO pQCD predictions and the QCDF predictions 
become narrow obviously because of the inclusion of the NLO contributions. 

For the five "QCD-Penguin" decays Bg rj^'^^'^ and KK decays, the enhancements 
due to the considered NLO contributions can be as large as (30 — 70)%, which are 
helpful to pin down the gap between the pQCD and the QCDF predictions. 

For the two "Electroweak- Penguin" decays Bg 7r^r]^'\ the NLO contributions are 
small. The pQCD predictions agree well with the QCDF predictions. 

For the "annihilation" decays Bg tt+tt^ and vr'^Tr*' decays, the NLO contributions 
are around 10% only. The pQCD predictions agree well with the measured value. 



For the considered thirteen Bg —>■ PP decays, only three of them, Bg 
K~^7T~ , K^K~ and 7r"'"7r~, have been measured experimentally with good preci- 
sion. It is easy to see that the consistency between the pQCD predictions for 
their branching ratios and the measured values will be improved effectively when 
the NLO contributions are included. 



B. CP- violating asymmetries 



Now we turn to the evaluations of the CP-violating asymmetries of Bg —>■ PP decays 
in pQCD approach. Restricting the final state f to have definite CP, the time-dependent 
decay width for the Bg ^ f decay can be written as [24 



/Art 



cosh 



+Hf sinh 



Acp cos(Am t) + Sf sin( Am t) 



(79) 



where Am = uih — mi > 0, F = {Vh + rL)/2) is the average decay widths, while AF = 
Th — Fi, is the difference of decay widths for the heavier and lighter mass eigenstates. 
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TABLE II: Direct CP asymmetries (in %) of Bs — > PP decays in the pQCD approach. The label 
LO means the leading order pQCD predictions, while +VC, +QL, +MP, as well as NLO means 
that the vertex corrections, the quark loops, the magnetic penguins, and all the above NLO 
corrections are added to the LO results, respectively. The errors in the table are defined in the 
context. For comparison, we also cite the leading-order pQCD predictions as given in Ref. 10l |. 
the QCDF results in Ref. jfi]. 



Mode 


Class 


LO 


+ vc 


+ QL 


+ MP 


NLO 


pQCD[10] 


QCDF[6] 






C 


65.0 


45.4 


38.6 


30.7 


OR 7+0.0+1.1+1.2 
'-0.1-2.0-1.5 


56.411° 


46.8ti| 




kS 


c 


-22.3 


-5.7 


-18.1 


-0.8 


„_ .4-2 (14-2 4-Wl ^ 
"J^-^-0.0-2. 5-0.3 


-19.9+53 


-36.6+22-3 






Pew 


0.3 


40.4 






/in /I +0.3+1.6+3.6 
^"-'•^-0.8-1.3-7.2 


_n 4+0-3 

"J-^-0.3 






ir^T]' 


Pew 


23.8 


52.5 






t-r, c;+2. 0+2.4+0.5 
^^•^-0.0-2.5-0.3 


20.6^1^ 


01 Q+27.2 

^' •°-28.8 




r]r] 


P 


-0.7 


-1.53 


1.2 


1.1 


n R+0.1+0.1+0.2 
"•"-0.0-0.0-0.0 


-0 6+°-^ 
"J-"3-0.5 


-1 6+2-4 


B^ ^ 


rjrj' 


P 


-L3 


-1.1 


-0.2 


-0.6 


n. 9+0.1 + 0.0+0.1 

"•^-0.1-0.0-0.1 


_i 3+0.1 
-'■•"-'-0.2 


n 4+0-5 
"•^-0.4 


B^ ^ 


Tj'r]' 


P 


1.9 


1.3 


1.0 


1.3 


1 /i+O.l+O.l+O.l 
^•^-0.1-0.1-0.2 


1 q+0.4 
^•^-0.5 


2 1+^-3 




K+TT- 


T 


25.7 


28.6 


25.7 


28.5 


nc; Q+4.1+1.5+2.7 
^'^•°-3. 8-0.7-5.0 


24 1+^-^ 
^^•-"--4.8 


c 7+15.6 
"^•'-15.3 


B^ ^ 




C 


66.9 


86.4 


-18.7 


-10.8 


00 n+3. 7+2.6+1. 6 
oo-U_4 5_6.7_i,6 


59.411,^5 


42+47 
^^-56 


B^ ^ 




ann 


-1.1 








p. 9+0.1+0.0+2.0 
'-'•^-0.0-0.0-1.5 






B^ 




ann 


-1.1 








"-'•^-1.5 


-1 2+^-2 




B^ ^ 


K+K- 


P 


-22.1 


-17.9 


-14.1 


-17.1 


1 r fj+l. 2+0.7+1. 3 
J- '-'•"-'-0.8-0.9-1.1 


-23.3tt[! 


^•"-11.6 


B^ ^ 




P 








0.3 





0.4 ±0.1 





0.9 ±0.4 



In the Bs system, we expect a much larger decay width difference: (Ar/r)^^ ~ — zu7o 
2l|. Besides A'^^^ the CP- violating asymmetry Sr and Hr can be defined as 



A dir 
J^CP 



with the parameter A 



|Ap-l 

i±iAr 



2Im[A] 



Hf 



2Re[A] 
1± lAP' 



(80) 



A = ?7/e 



A{Bs^f)' 



(81) 



where r]f is ±1(— 1) for a CP-even(CP-odd) final state f and e = arg[— Vj^V^^] is very small 
in size. 

If we neglect the very small parameter e, the CP-violating asymmetries can be written 
explicitly as 

jdir ^ 22: sin a sin 6 

1 ± 2zcosacos(5 + ' 



Sf 
Hf 



sin(27) ± z"^ sin(27 + 2a) + 2z cos 5 sin(a ± 27) 
1 ± ^2 -|- 22; cos 6 cos a ' 
2z cos((5) cos(a ± 27) ± cos(27) ± z"^ cos(2q; ± 27) 
1 + z"^ + 2z cos 6 cos a 



m 
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for the decays relevant to the h ^ d transition, and 

2z' sin 7 sin 5 



A dir 



S 



f 



1 + 2z' cos 7 cos 6 + z'"^ 
sin(27) + 2z' cos 6 sin 7 
1 + _)_ 2z' cos 5 cos 7 ■ 



^ _ z'^ + 22;' cos 5 cos 7 + cos (27) 
^ ~ l + z'2 + 22' cos 5 cos 7 ' ^ 

for the case of 6 s transition. 

The pQCD predictions for the direct CP asymmetries A'^^^, the mixing-induced CP 
asymmetries Sf and Hf of the considered PP decays are hsted in Table [Til and 

Table IIIII In both tables, the label LO means the leading order pQCD predictions, and 
the labels +VC, +QL, +MP, as well as NLO mean that the vertex corrections, the quark 
loops, the magnetic penguins, and all the above NLO corrections are added to LO results, 
respectively. As a comparison, the LO pQCD predictions as given in Ref. [lol] are listed 
in Table [IT] and Table IIIII In Table Ull the QCDF predictions for direct CP-violating 
asymmetries as given in Ref. are also shown. The corresponding errors of the previous 
LO pQCD predictions and QCDF predictions are the combined errors: the individual 
errors as given in Refs. 0, [lO] are added in quadrature. The errors of our NLO pQCD 
predictions for CP-violating asymmetries are defined in the same way as those for the 
branching ratios. 

For the experimental measurements, there is only one measured CP asymmetry as 



reported by CDF Collaboration |22|: 

A'^jpiBs K+n-) = 0.39 ± 0.17. (84) 

But more data will become available soon when the LHC starts its physics running. 

From the pQCD predictions and currently available experimental measurements for 
the CP violating asymmetries of the thirteen B PP decays, one can see the following 
points: 

• The LO pQCD predictions obtained in this paper agree very well with those as 
given in Ref. [13] • For Bg K^rj and tt^t] decays, the LO pQCD predictions can be 
changed significantly by the inclusion of the NLO contributions. For other decays, 
the NLO contributions are small or moderate in size. The pQCD predictions are 
in general consistent with those in QCDF approach, but much larger than the later 
one for Bs K^, rc^rj' and K^tc^ decays. 

• For the "Tree" dominated decay Bg K^tt^ , the pQCD prediction for the direct 
CP asymmetry is A'^U^s ^ K+tt-) = 0.30 ±0.06, which agrees very well with the 
experimental measurement as given in Eq. flH^ . The QCDF prediction, however, is 
about —0.07 ± 0.16 and much different from the measured value. 

• For the four "QCD-penguin" decays Bg r]^''>r]^'^ and K^K^ decays, analogous to 
the QCDF predictions, the LO and NLO pQCD predictions for both A^p and 5*/ 
are all very small in size. 
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TABLE III: The mixing-induced CP asymmetries (in %) Sj and Hf (the second row). The label 
LO means the leading order pQCD predictions, and the lables +VC, -|-QL,-|-MP, as well as NLO 
mean that the vertex corrections, the quark loops, the magnetic penguins, and all the above 
NLO corrections are added to LO results, respectively. The errors of the entries are defined in 
the context. As a comparison, the LO pQCD predictions as given in Ref. [lO(] are also listed. 



Mode 


Class 


LO 


+VC +QL +MP 


NLO 


pQCD[10] 






C 


-37 
-67 


-89 
-57 


-92 
6 


-90 
32 


1 e+O+7+3 
-i8_2-ll-5 
18+2+18+15 

-■-"-o-s-g 


-43li 






C 


-67 


-59 


-44 


-53 


-46+1+12+1 

^"-0-23-0 


-68^^ 








-70 


-80 


-88 


-85 


00+0+20+1 
°^-l-6-0 


-70+? 




TT^r] 


Pew 


18 


28 


— 


— 


90+2+3+4 
^^-3-1-4 


17^1^ 








98 


87 


— 


— 


Q7+1+1+4 
0'-l-l-2 


99i^ 


B's^ 


ttS 


Pew 


-25 


-18 


— 


— 


10+1+12+1 
■'-°-0-23-0 


-17+8 
I'-g 








94 


83 


— 


— 


oQ+3+17+1 


96^2 


B's^ 


rjiTj 


P 


3 





1 


1 


9+0+0+0 
^-0-0-0 




B's^ 


rjT]' 


P 


100 
4 


100 
3 


100 
4 


100 
3 


100+°+°+° 

/I +0+0+0 
^-0-0-0 


100+° 

4t° 




rj'r]' 


P 


100 

4 


100 
6 


100 
6 


100 

6 


100+°+°+° 

r-+0+0+0 

^-1-1-0 


100+° 

At\ 








100 


100 


100 


100 


- ^^4-04-04-0 

iUU_Q_Q_Q 


lOOlg 






C 


-55 
-50 


-25 
-44 


-98 
-8 


-97 
-20 


1+8+4+3 
^J^-9-8-5 
90+0+19+3 
^•5-1-18-7 


_fi1+24 

-52l?3 




K+K- 


P 


24 


20 


22 


20 


99+2+2+2 
^^-2-1-2 










95 


96 


96 


96 


Qf;+4+0+0 
yO-3-1-0 


93^3 


B's^ 




P 






0.4 




n /I +0+0+0 
•-"•^-o-o-o 


4 












100 




100+°+°+° 


100 


B's^ 
B's^ 




ann 
ann 


9.5 
99.5 
9.5 








^-0-0-0 

1 nn+0+0+0 

iUU_o_o_o 
Q 1+0.1+0.5+0.3 
"■1-0.3-0.7-0.0 


99+° 








99.5 








100+°+°+° 

-LUU_Q_Q_Q 


99+° 



• For the "annihilation" decays Bs — > tt+tt^ and 7r°7r°, the pQCD predictions for the 
direct CP-violating asymmetries are very small in size, while Sf is around 10% and 
Hf ~ 1. 



VI. SUMMARY 

In this paper, we calculated the partial NLO contributions to the branching ratios and 
CP-violating asymmetries of 5° — »• PP decays. Here the NLO contributions from the 
QCD vertex corrections, the quark-loops and the chromo- magnetic penguins are included. 

From our calculations and phenomenological analysis, we found the following results: 
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The LO pQCD predictions for the branching ratios and CP-violating asymmetries 
of Bs — >■ PP decays as presented in Ref. [lO| are confirmed by our independent 
calculation. 



• For branching ratios, the effects of the considered NLO contributions are varying 
from small to significant for different decay mode. For the three measured decays 
Bg K^n^ , K^K^ and tt+tt^, for example, the consistency between the pQCD 
predictions and the measured values are improved effectively due to the inclusion 
of the considered NLO contributions. For the three "Color-suppressed" decays, for 
instance, the NLO enhancement can be significant, from ~ 50% for the Bg — >■ K'^t] 
and X°7r° decays to ~ 170% for the Bg — > K'^r]' decay, to be tested by forthcoming 
LHC experiments. 

• As for the CP-violating asymmetries, the LO pQCD predictions for Bg — » K^t] and 
7r°?7 decays could be changed significantly by the inclusion of the NLO contributions. 
For other decays, the NLO contributions are small or moderate in size. For Bg 
K~^n~ decay, the pQCD prediction for the direct CP asymmetry is A^plBg 
K^Ti~) = 0.26 ± 0.06, which agrees very well with the measured value AQp{Bg 



• In this paper, only the partial NLO contributions in the pQCD approach have been 
taken into account. The still missing pieces relevant with the emission diagrams, 
hard-spectator and annihilation diagrams should be evaluated as soon as possible. 
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APPENDIX A: RELATED FUNCTIONS 

We show here the hard function hi and the Sudakov factors Sab,cd,ef,ghif) appeared in 
the expressions of the decay amplitudes in Sec. [ITIlandllVl The hard functions hj) 
are obtained by making the Fourier transformations of the hard kernel H^^\ 



K+n-) = 0.39 ±0.17. 



(Al) 



hn{xi,bi,b2) = [6'(&2 - &i)Ko(MB^v/i7a^62)Io(MB^V'xIx^6i) 



+^(61 - 62)Ko(MB^y^6i)Io(MB,v/^62)] 




(A2) 
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/ina(a;i, 61, &2) = Oibi - 62)Ko (^2^/0:2(1 - X3)MbM^ Io (^2^2(1 - a;3)Mij>2) 
+e{b2 - 6i)Ko (Va;2(l-a;3)MB,&2) Iq {^^/ X2{1 - X3)MbM 
xKo - (1 - a;i - X2)x3MbM 



(A4) 



^na(a;i, &1,&2) 



0(61 - 62)Ko (zVa;2(l-a;3)MB,6i) Io {l^J X2{1 - X3)MB^b 
+e{b2 - &i)Ko (l^Jx2{l-X3)MBb^ Io fzVx2(l-a;3)MB,6i 



X 



Mb^ a/ (x2 - xi)(l - X3)6ij , for xi- X2 < 0, 
Ko (^Mb, a/(x2 - xi)(l - X3)6ij , for X2- xi> 0, 



(A5) 



bi) 



7r/2 



dOtanO ■ Jo (A/a^MB^6itan6') Jo i^/x^MB,b2tan9) 
■Jo (v/a^MB^63tan6') , 



(A6) 



h'Jxi,bi) 



where 



/•'r/2 

-5t(xi)iro (ViI^MB,&3) • / detank- Jo (v^Mb, 61 tan ^) 
■Jo i\/xiMB,b2 tan 61) Jo (v/a^MB^63 tan 6) 



^ ^ f [Jo (Va;2 - xiMbM) + ^A'o (V2;2 - xiMbM)] , xi < X2, ^j^^^ 



Ko {y/Xi - X2MbM) 



Xi > X2, 



ITT 



— [iYo{x) + Jo{x)] 



(A8) 



with Ko, Io and Jo are the Bessel functions j25|]. And the threshold resummation form 
factor St{ be found in Ref. 261. 

The Sudakov factors appeared in Eqs. (13711461) are defined as 



Sab{t) = s { Xi^^, bi] + s{ Xs^^, 63 ) + s ( 0:3^^, 63 



5 

+3 



1/61 



/i 



(A9) 
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Scd{t) = S [ ^l—J^^ ^ij + ^ ^^j + S (^^2^^, &2j + S hi 



+ S 1 X3-^, b,j + - 



1/fel il/fe2 



(AlO) 



Sef{t) = S { X2—^, 02] + S { X2—^, 02] + S { X3—^, 63 



V2 



V2' 



V2 



1/62 



2 r rf/.^^^^, (All) 



+ S X3 



V2 

rriB 



V2 



S I Xi—^, h] + S { X2—^, O2 1 + S I X2—^, 02] + S { X3—^, 62 



V2 



V2 



V2 / 3 Ji/h^ /i Ji/bj 



(A12) 



where the quark anomalous dimension 7^ = —ag/iT and the function s{Q,b) is given as 
13, 



siQ,b) 



dfi 

1/6 



ln(^)A(a(^)) + 5(a,(/i)) 



(A13) 



with 



A 



3 TT 

2 a, 



+ 



7B 



5 = ^In 

3 TT 



67 TT^ 10 2^ / e 

A^/- + -/?oln — 

9 3 27 ^ 3^° V 2 



TT 



(A14) 



where 7^; = 0.57722 ■ ■ ■ is the Euler constant, an Nf is the number of active quark flavors. 
The hard scales appeared in the above equations take the form of 



max{y/x^MB,, l/bi, I/63} , 
max{y/x^MB,, l/bi, I/&3}, 

max Ia/xiXsMb^, ^y\ 1 - xi - X2 \ x^Ms^.l/bi, I/62I 



max j^XiXsM^^, xi - X2 | a;3Mij^, I/61, 1/62 

max {a/I^Mb^, 1/62, 1/63} , 
max{v/a^MB^, I/62, 1/&3} , 

max |a/x^Mb^, ^/l - (1 - xi - X2)x3Mb,, I/61, 1/&2I 
max |a/x2X3Mb^, a/| xi - X2 | X3Mij^, I/61, 1/62 



(A15) 
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They are chosen as the maximum energy scale appearing in each diagram to kill the large 
logarithmic radiative corrections. 
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